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Abstract 

We present the Flavour Les Houches Accord (FLHA) which specifies a set of con- 
ventions for flavour-related parameters and observables. The FLHA uses the generic 
SUSY Les Houches Accord (SLHA) file structure. It defines the relevant Standard 
Model masses, Wilson coefficients, decay constants, bag parameters, flavour observ- 
ables, etc. The accord provides a universal and model-independent interface between 
codes evaluating and/or using flavour-related observables. 
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1 Introduction 



In addition to the increasing number of refined approaches in the literature for calculat- 
ing flavour-related observables, advanced programs dedicated to the calculation of such 
quantities, e.g. Wilson coefficients, branching ratios, mixing amplitudes, renormalisation 
group equation (RGE) running including flavour effects have recently been developed [IHI]. 
Flavour-related observables are also implemented by many other non-dedicated public codes 
to provide additional checks for the models under investigation [5TTTB] . The results are often 
subsequently used by other codes, e.g. as constraints on the parameter space of the model 
under consideration [T1HTT] . 

At present, a small number of specialised interfaces between the various codes exist. 
Such tailor-made interfaces are not easily generalised and are time-consuming to construct 
and test for each specific implementation. A universal interface would clearly be an ad- 
vantage here. Some time ago a similar problem arose in the context of Supersymmetry 
(SUSY). The solution took the form of the SUSY Les Houches Accord (SLHA) [THEE], 
which is nowadays frequently used to exchange information between SUSY related codes, 
such as values for the soft SUSY-breaking parameters, particle masses and mixings, branch- 
ing ratios etc. The SLHA has proved to be a robust solution, allowing information to be 
exchanged between different codes via ASCII files. The detailed structure of these input 
and output files is described in Refs. [T5}[T9"]. While the first definitions, SLHAl [T8], con- 
cerned the simplest version of the MSSM, the SLHA2 [19] also included definitions for the 
case of complex parameters, i?-parity violation, non-minimal flavour violation etc. 

The goal of this article is to exploit the existing organisational structure of the SLHA 
and use it to define an accord for the exchange of flavour related quantities, which we 
refer to as the "Flavour Les Houches Accord" (FLHA). In brief, the purpose of this Ac- 
cord is thus to present a set of generic definitions for an input/output file structure which 
provides a universal framework for interfacing flavour-related programs. Furthermore, the 
standardised format will provide the users with a clear and well-structured result that could 
eventually be used for other purposes. We stress that in cases of ambiguity in the precise 
definition of a given quantity, it is the responsibility of the authors of the specific program 
to document in detail the definition they use. 

The structure is set up in such a way that the SLHA and the FLHA can be used together 
or independently. Obviously, some of the SLHA entries, such as measured parameters in the 
Standard Model (SM) and the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements, are 
also needed for flavour observable calculations. Therefore, a FLHA file can indeed contain 
a SLHA block if necessary. For this reason, and also for the sake of clarity, the new FLHA 
block names start with "F" . Also, in order to avoid any confusion, the SLHA blocks are not 
modified or redefined in the FLHA. Instead, if a block needs to be extended to meet the 
requirements of flavour physics, a new "F" block is defined. 

Note that different codes may technically achieve the FLHA input /output in different 
ways. The details of how to "switch on" the FLHA input/output for a particular program 
should be described in the manual of that program and are not covered here. For the 
SLHA, libraries have been developed to permit an easy implementation of the input/output 
routines [20]. In principle these programs could be extended to include the FLHA as well. 
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It should be noted that, while the SLHA was developed especially for the case of SUSY, 
the FLHA is, at least in principle, model independent. While it is possible to indicate the 
choice of model in a specific block, the general structure used for the information exchange 
can be applied to any model. 

2 Definitions of the interfaces 

The FLHA input and output files are described in this section. 
2.1 General structure 

Following the general structure for the SLHA [18j[l9] we assume the following: 

• All quantities with dimensions of energy (mass) are in GeV (GeV/c 2 ). 

• Particles are identified by their PDG particle codes. See Appendix [A] for lists of these, 
relevant for flavour observables. 

• The first character of every line is reserved for control and comment statements. The 
first character of data lines should be empty. 

• In general, a formatted output should be used for write-out, to avoid "messy-looking" 
files, while a free format should be used for read-in, to avoid misalignment etc. leading 
to program crashes. 

• Read-in should be performed in a case-insensitive way, again to increase stability. 

• The general format for all real numbers is the FORTRAN format E16.s(E The large 
number of digits is used to avoid any possible numerical precision issue, and since 
it is no more difficult for e.g. the spectrum calculator to write out such a number 
than a shorter version. For typed input, this merely means that at least 16 spaces 
are reserved for the number, but e.g. the number 123.456 may be typed in "as is". 
See also the example file in Appendix [G] 

• A "#" mark anywhere means that the rest of the line is intended as a comment and 
should be ignored by the reading program. 

• Any input and output is divided into sections in the form of "blocks" . 

• To clearly identify the blocks of the FLHA, the first letter of the name of a block 
is an "F". There are two exceptions to this rule: blocks borrowed from the SLHA, 
which keep their original name, and blocks containing imaginary parts, which start 
with "IMF". 

^^EIG.S: a 16-character wide real number in scientific notation, whereof 8 digits are decimals, e.g., 
"-0.12345678E+000". 
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• A "BLOCK Fxxxx" (with the "B" being the first character on the line) marks the 
beginning of entries belonging to the block named "Fxxxx". For instance, "BLOCK 
FMASS" marks that all following lines until the next "BLOCK" statement contain mass 
values, to be read in a specific format, intrinsic to the FMASS block. 

• The FLHA is designed to be compatible with SLHA (and with any future accords that 
follow the same general structure). For generic models, e.g. ones for which an accord 
does not yet exist, we provide a fall-back solution by generalising the SLHA block 
MODSEL. This takes the form of a new FLHA-specific block, FMODSEL, with entries as 
described below. For models for which a more specific accord does exist, e.g. SLHA, a 
complete model of that type can be specified instead, such that FMODSEL is redundant 
and should in most cases be absent. In case it is present and there are conflicts, the 
most specific accord takes precedence, such that internal consistency of the accord 
with the smallest scope is guaranteed. In the case of MODSEL vs. FMODSEL, MODSEL 
would thus take precedence, however it would be prudent to issue warnings if both 
are present, especially if their contents differ. Also note that although there is the 
possibility to reuse more specific accords, one is not forced to do so. Thus, even 
for a SUSY model, an FLHA user who does not wish to deal with an entire SLHA 
spectrum could just use FMODSEL. In that case, the file would be treated just like any 
other FLHA file, i.e. without any reference to the SUSY-specific parts of the SLHA. 
The only difference from a practical perspective would be that SLHA-specific tools 
would then not be able to process the file correctly. 

• The order of the blocks is arbitrary, although it is in general good practice to put 
FMODSEL or MODSEL near the beginning. 

Further definitions can be found in section 3 of Ref . [18] . 

The following general structure for the FLHA file is proposed: 

• BLOCK FCINFO: Information about the flavour code used to prepare the FLHA file. 

• BLOCK FMODSEL: Basic information about the underlying model used for the calcula- 
tions, for generic models. This is the only place where "model dependent" information 
can be found. In the case of SUSY models with complete SLHA spectra, the SLHA 
BLOCK MODSEL is used instead and overrides FMODSEL if both are present. 

• BLOCK SMINPUTS: Measured values of SM parameters used for the calculations. 

• BLOCK VCKMIN: Input parameters of the CKM matrix in the Wolfenstein parameteri- 
sation. 

• BLOCK UPMNSIN: Input parameters of the PMNS neutrino mixing matrix in the PDG 
parameterisation. 

• BLOCK VCKM: Real part of the CKM matrix elements. 
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• BLOCK IMVCKM: Imaginary part of the CKM matrix elements. 

• BLOCK UPMNS: Real part of the PMNS matrix elements. 

• BLOCK IMUPMNS: Imaginary part of the PMNS matrix elements. 

• BLOCK FMASS: Masses of quarks, mesons, hadrons, etc. 

• BLOCK FLIFE: Lifetime (in seconds) of mesons, hadrons, etc. 

• BLOCK FCONST: Decay constants. 

• BLOCK FCONSTRATIO: Ratios of decay constants. 

• BLOCK FBAG: Bag parameters. 

• BLOCK FWCOEF: Real part of the Wilson coefficients. 

• BLOCK IMFWCOEF: Imaginary part of the Wilson coefficients. 

• BLOCK FOBS: Prediction of flavour observables. 

• BLOCK FOBSERR: Theory error on the prediction of flavour observables. 

• BLOCK FOBSSM: SM prediction for flavour observables. 

• BLOCK FDIPOLE: Prediction of electric and magnetic dipole moments. 

• BLOCK FPARAM: Process dependent variables and parameters. 
More details on each block are given in the following. 

2.2 Definition of the blocks 

The FLHA input and output blocks are described in the following. 
BLOCK FCINFO 

Flavour code information, including the name and the version of the program: 

1 : Name of the flavour calculator 

2 : Version number of the flavour calculator 

Optional warning or error messages can also be specified: 

3 : If this entry is present, it means warnings were produced by 

the flavour calculator. The resulting file may still be used. 
The entry should contain a description of the problem (string). 

4 : If this entry is present, it means errors were produced by the 

flavour calculator. The resulting file should not be used. The 
entry should contain a description of the problem (string). 

A string is as usual a sequence of ASCII characters. This block is purely informative, and 

is similar to BLOCK SPINFO in the SLHA. 
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BLOCK MODSEL 



This block is part of the SLHA. Its presence in an FLHA file signals that the file contains, 
in addition to the FLHA-specific blocks, a complete and self-consistent SLHA spectrum, 
for which MODSEL defines the switches and options used. It is defined as in the SLHA2. For 
non-SLHA models, see instead BLOCK FMODSEL below. Note that, if the user does not wish 
to provide a complete set of SLHA blocks, MODSEL should not be used. Instead, see BLOCK 
FMODSEL below. 



BLOCK FMODSEL 



In the case of non-SLHA models this block provides switches and options for model selection. 
It is similar to the SLHA2 BLOCK MODSEL, but extended to allow more flexibility. It is 
advised that one only uses BLOCK MODSEL for SLHA models and BLOCK FMODSEL for non- 
SLHA models, to avoid double definitions. 

1 : Choice of model. By default, a minimal type of model will 
always be assumed. Possible values are (note: when giving a 
complete SLHA spectrum, use BLOCK MODSEL instead. The 
options here are only intended to cover the cases when a 
complete SLHA spectrum is not provided.): 

SM 

General MSSM simulation 
(m)SUGRA model 
(m)GMSB model 
(m)AMSB model 



-1 


1 
2 
3 
4 
30 
31 
32 
33 
34 
35 
99 



General THDM 
THDM type I 
THDM type II 
THDM type III 
THDM type IV 



other model. 
99 



This choice requires a string given in the entry 



(Default=0) CP violation. Switches defined are: 

: CP violation is completely neglected. No information on the 

CKM phase is used. 

1 : CP is violated, but only by the standard CKM phase. All 

other phases are assumed zero. 

2 : CP is violated. Completely general CP phases allowed. 
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6 : (Default=0) Flavour violation. Switches defined are: 

: No flavour violation. 

1 : Quark flavour is violated. 

2 : Lepton flavour is violated. 

3 : Lepton and quark flavour is violated. 

99 : a string that defines other models is used only if entry 1 is 
given as 99, otherwise it is ignored. 

The definition of the different Two-Higgs Doublet Model (THDM) types is given in 
Appendix IB1 

Private blocks can also be constructed by the user, for instance BLOCK MYMODEL, to 
contain parameters specific to other models. It is advised that all such relevant model 
parameters are provided in this way. Due to the user-defined specific structure of these 
blocks it is not required that they are universally recognised. 



BLOCK SMINPUTS 

In general, the spectrum of the SM particles plays a crucial role in flavour physics. Conse- 
quently, experimental measurements of masses and coupling constants at the electroweak 
scale are required. The block containing these quantities in the SLHA is SMINPUTS. We 
borrow this block from SLHA as it is, and reproduce it here for completeness. 

It is also important to note that experimental results for all quantities available at 
present, e.g. a s or the running bottom quark mass, are clearly obtained based on the 
assumption that the SM is the underlying theory. Extending the field content of the SM to 
that of a New Physics (NP) Model means that the same measured results would be obtained 
for different values of these quantities. However, since these values are not known, all 
parameters contained in the block SMINPUTS should be the "ordinary" ones obtained from 
SM fits, i.e. with no NP corrections included. Any flavour code itself is then assumed to 
convert these parameters into ones appropriate to an NP framework. 

It should be noted that some programs have hard-coded defaults for some of these pa- 
rameters, hence only a subset may sometimes be available as free inputs. The parameters 
are: 

1 : a~^(mz) MS , inverse electromagnetic coupling at the Z pole in 

the MS scheme (with 5 active flavours). 

2 : Gf-i Fermi constant (in units of GeV -2 ). 

3 : a s (mz) MS , strong coupling at the Z pole in the MS scheme 

(with 5 active flavours). 

4 : mzi pole mass. 

5 : rrib(mb) MS , bottom quark running mass in the MS scheme 

(with 5 active flavours). 

6 : m t , top quark pole mass. 

7 : m T , tau pole mass. 
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BLOCK VCKMIN 



This block is strictly identical to the SLHA2 BLOCK VCKMIN. The parameters are: 



1 


: A 


2 


: A 


3 


■ P 


4 


: fj 



We use the PDG definition, Eq. (11.4) of Ref. [21], which is exact to all orders in A, and 
also the PDG parameterisation for the phase convention. For the output we use the SLHA2 
blocks VCKM and IMVCKM which do not rely on any specific convention. 

BLOCK UPMNSIN 

This block is strictly identical to the SLHA2 BLOCK UPMNSIN. The parameters are: 



1 


: 012 


2 


'■ 023 


3 


'■ &13 


4 


: 5 


5 


: a,\ 


6 


: a 2 



We use the PDG parameterisation, Eq. (13.30) of Ref. [21j . All the angles and phases should 
be given in radians. For the output we use BLOCK UPMNS and BLOCK IMUPMNS, exactly in 
the same way as in SLHA2. 

BLOCK FMASS 

The block BLOCK FMASS contains the mass spectrum for the particles involved, in addition to 
the SLHA BLOCK MASS which only contains pole masses and to the SLHA BLOCK SMINPUTS 
which contains quark masses. If a mass is given in two blocks the block FMASS overrules the 
other blocks. In FMASS we specify additional information concerning the renormalisation 
scheme as well as the scale at which the masses are given and thus allow for larger flexibility. 
The standard for each line in the block should correspond to the FORTRAN format, 

(1x,I9,3x,1P,E16.8,0P,3x,I2,3x,1P,E16.8,0P,3x, '#' ,lx,A), 

where the first nine-digit integer should be the PDG code of a particle, followed by a double 
precision number for its mass. The next integer corresponds to the renormalisation scheme, 
and finally the last double precision number indicates the energy scale (0 if not relevant). 
An additional comment could be given after #. 
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The schemes are defined as follows: 






pole 


1 


MS 


2 


DR 


3 


IS 


4 


kin 


5 





For the definition of DR scheme see [221123] . 
BLOCK FLIFE 

The block BLOCK FLIFE contains the lifetimes, in seconds, of mesons and hadrons. The 
standard for each line in the block should correspond to the FORTRAN format 

(lx,I9,3x,lP,E16.8,0P,3x, '#' ,lx,A), 

where the first nine-digit integer should be the PDG code of a particle and the double 
precision number its lifetime. 

BLOCK FCONST 

The block BLOCK FCONST contains the decay constants in GeV. The renormalisation scheme 
and scale are also specified to allow for large flexibility. The standard for each line in the 
block should correspond to the FORTRAN format 

( lx , 19 , 3x , 12 , 3x , IP , E16 . 8 , OP , 3x , 12 , 3x , IP , E16 . 8 , OP , 3x , ' # ' , lx , A) , 

where the first nine-digit integer should be the PDG code of a particle, the second integer 
the number associated with the decay constant, the double precision number the value of 
this decay constant, the following integer stands for the renormalisation scheme (see below) 
and finally the last double precision number indicates the renormalisation scale Q (0 in 
the case of renormalisation group invariant parameters). An additional comment could be 
given after #. 

The decay constants for the mesons that are used most often, which have several decay 
constants associated with them, are defined as: 

321 : K + 

1 : f K in GeV 
11 : h K in GeV 3 
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221 



1 
2 
11 
12 

213 : p(770) 4 
1 

11 

223 : w(782) 
1 

2 

11 

12 



/' in GeV 
/• in GeV 
h\ in GeV 3 
h 8 in GeV 3 



f P m GeV 
fj in GeV 

P P in GeV 
fp in GeV 
fj>* in GeV 
/J> s in GeV 



More details, and definitions for the decay constants (/, h, etc.) can be found in Ap- 
pendix [D] 

Finally, the schemes are defined as follows: 

: renormalisation group invariant (with renormalisation scale 
set to 0) 



1 

2 
3 
4 



MS 

MOM 

SMOM 



BLOCK FC0NSTRATI0 

The block BLOCK FC0NSTRATI0 contains the ratios of decay constants, which often have 
less uncertainty than the decay constants themselves. The ratios are specified by the two 
PDG codes in the form f(codel)/f(code2). The standard for each line in the block should 
correspond to the FORTRAN format 

(1x,I9,3x,I9,3x,I2,3x,I2,3x,1P,E16.8,0P,3x,I2,3x,1P,E16.8,0P,3x, '#' ,lx,A), 

where the two nine-digit integers should be the two PDG codes of particles, the third 
and fourth integers the numbers of the decay constants, which correspond to the second 
index of the entry in BLOCK FCONST, the double precision number the ratio of the decay 
constants and the following integer stands for the renormalisation scheme and finally the 
last double precision number indicates the renormalisation scale Q in the same way as in 
BLOCK FCONST. 
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BLOCK FBAG 



The block BLOCK FBAG contains the bag parameters. The renormalisation scheme as well 
as the scale at which the bag parameters are evaluated are also specified to allow for large 
flexibility. The standard for each line in the block should correspond to the FORTRAN 
format 

( lx , 19 , 3x , 12 , 3x , IP , E16 . 8 , OP , 3x , 12 , 3x , IP , E16 . 8 , OP , 3x , ' # ' , lx , A) , 

where the first nine-digit integer should be the PDG code of a particle, the second integer 
the number associated with the bag parameter, the double precision number the value of 
the bag parameter, the following integer stands for the renormalisation scheme and finally 
the last double precision number indicates the renormalisation scale Q (0 in the case of 
renormalisation group invariant bag parameters). An additional comment could be given 
after #. 

Bag parameters for B® — B d (and similarly for B° s — B®, K° — K° and D° — D°) mixing 
are listed as: 

511 : B° d 

1 : B 1 (Q) 

2 : B 2 (Q) 

3 : B 3 (Q) 

4 : B 4 {Q) 

5 : B 5 (Q) 

The schemes are defined in the same way as in BLOCK FCONST and the Bag parameter 
definitions and normalisations can be found in Appendix [El 



BLOCK FWCOEF Q= . . . 

The block BLOCK FWCOEF = • • • contains the real part of the Wilson coefficients at the 
scale Q. The Wilson coefficients are calculated for the required set of operators. A list of the 
most relevant effective operators is given in Appendix O and we recommend the user to use 
the normalisation as given in that appendix (the detailed information about the scheme and 
the normalisation of the effective interactions has to be provided in the program manual, 
and a short comment can be given in the FLHA file). 

Note that there can be several such blocks for different scales Q. 

The coefficients C^j at order k in the perturbative expansion in a(fi) (j = e) or a s (fi) 
(J = s), have to be specified individually. The following convention for the perturbative 
expansion is used: 



^c^ ) + a M^lc^) + .... (i) 
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name 


PDG code 


two-digit number 


name 


PDG code 


two-digit number 


d 


1 


01 


e 


11 


11 


u 


2 


02 




12 


12 


s 


3 


03 




13 


13 


c 


4 


04 


v„ 


14 


14 


b 


5 


05 


T 


15 


15 


t 


6 


06 


V T 


16 


16 


Eg Q 




07 






17 






08 






18 



Table 1: PDG codes and two-digit number identifiers of quarks and leptons. The summa- 
tions are over active fermions. 



Therefore the couplings should not be included in the Wilson coefficients. 

The first two entries on each line in BLOCK FWC0EF should consist of two integers defining 
the fermion structure of the operator and the operator structure itself. These two numbers 
do not describe the operator fully, including normalisation etc. , but together act as a 
unique identifier for any possible Wilson coefficient. Consequently, the user has to take 
care that a consistent normalisation including prefactors etc. is used in specifying the 
Wilson coefficients. The most relevant operators are listed in Appendix As an example, 
for the operator 1} 

1 = (7s^T a P L c)(crT a P L b) , (2) 

the definition of the two numbers is given as follows. The fermions are encoded by their 
PDG code in two-digit form, which is the same for particles and antiparticles, as given in 
Table [TJ Correspondingly, the first integer number defining Oi, containing the fermions 
sccb, is given by 03040405. The various operator sub-structures are defined in Table |2j 
Correspondingly, the second integer number defining Ox, containing the operator structure 
^^T a P L ^T a P L is given by 6161. 

A few more rules are needed for an unambiguous definition. 

• If part of an operator appears without fermions (as it is possible, e.g., for F^ w ) it 
should appear right-most, so that the encoded fermions correspond to the left-most 
part of the operator. 

• In the case of a possible ambiguity, for instance 0\ = {s^ ll T a PLc){c r j^ l T a PiJj) cor- 
responding to 03040405 6161 and 0\ = (c r y^T a P L b)(s^T a P L c) corresponding to 
04050304 6161, the "smaller" number, i.e. in this case 03040405 6161, should be 
used. 

The third entry specifies the order of the perturbative expansion in Eq. (JTJ). The infor- 
mation about the order is given by a two-digit number xy, where x indicates 0(a x ) and y 
indicates 0{a y s ), and indicates Cf \ e.g.: 

00 : Cf (Q) 
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operator 


number 


operator 


number 


operator 


number 


1 


30 




50 




70 


Pl 


31 


P L T a 


51 


P L Sij 


71 


Pr 


32 


P R T a 


52 


PROij 


72 


r 


33 


^T a 


53 




73 


75 


34 




54 


75% 


74 




35 




55 




75 




36 




56 




76 




37 




57 




77 


YPl 


41 


YT a P L 


61 


-fSijPi, 


81 


YPr 


42 


l*T a P R 


62 


rkjPR 


82 




43 


a^T a P L 


63 


O^kjPL 


83 


O^Pr 


44 


a^T a P R 


64 


O^kjPR 


84 


^YYPl 


45 


^^Y^ p T a P L 


65 


n u Y^p L 


85 




46 




66 


^rrkjPR 


86 


F 


22 


" iiv 


21 







Table 2: Two-digit number definitions for the operators. T a (a = 1 . . . 8) denote the SU(3)c 
generators, P LiR = ±(1 =F7 5 ), and (T a ) i:j (T a ) k i = \(5 u 5 kj - l/N c 5ijS k i), where i,j,k,l are 
colour indices. 



01 


: Ct\Q) 


02 


■ C%(Q) 


10 


■ C%{Q) 


11 


■ C®(Q) 


99 


: total 



The Wilson coefficients can be provided either via separate New Physics and SM contri- 
butions, or as a total contribution of both New Physics and SM, depending on the code 
generating them. To avoid any confusion, the fourth entry must specify whether the given 
Wilson coefficients correspond to the SM contributions, New Physics contributions or to 
the sum of them, using the following definitions: 

: SM 

1 : NP 

2 : SM+NP 

The New Physics model is the model specified in the BLOCK FM0DSEL. 

The standard for each line in the block should thus correspond to the FORTRAN format 

(lx,I8,lx,I4,3x,I2,3x,Il,3x,lP,E16.8,0P,3x, '#' ,lx,A), 

where the eight-digit integer specifies the fermion content, the four-digit integer the operator 
structure, the two-digit integer the order at which the Wilson coefficients are calculated 
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followed by the one-digit integer specifying the model, and finally the double precision 
number gives the real part of the Wilson coefficient. 



BLOCK IMFWCOEF Q= . . . 

The block BLOCK IMFWCOEF contains the imaginary part of the Wilson coefficients at the 
scale Q. The structure is exactly the same as that of the BLOCK FWCOEF. 



BLOCK FOBS 

The block BLOCK FOBS contains the predictions for the flavour observables. The structure 
of this block is based on the decay table in SLHA format. The decay is defined by the 
PDG number of the parent, the type of observable, the value of the observable, the energy 
dependence q for the observable (0 if not relevant), number of daughters and PDG IDs of 
the daughters. 

The types of the observables are defined as follows: 

1 : Branching ratio 

2 : Ratio of the branching ratio to the SM value 

3 : Asymmetry - CP 

4 : Asymmetry - isospin 

5 : Asymmetry - forward-backward 

6 : Asymmetry - lepton-flavour 

7 : Mixing 

8 : ... 

The standard for each line in the block should correspond to the FORTRAN format 



(1x,I9,3x,I2,3x,1P,E16.8,0P,3x,1P,E16.8,0P,3x,I1,3x,I9,3x,I9,3x, . . . ,3x, , lx 

where the first nine-digit integer should be the PDG code of the parent decaying particle, 
the second integer the type of the observable, the double precision number the value of 
the observable, the next double precision number the energy scale q, the next integer the 
number of daughters, and the following nine-digit integers the PDG codes of the daughters. 
It is strongly advised to give the descriptive name of the observable as a comment. More 
details and definitions regarding the meson mixing are provided in Appendix [F] 

For more specific conventions for each observable, the user is encouraged to refer to the 
manual of the flavour calculator. 



BLOCK FOBSERR 

The block BLOCK FOBSERR contains the theoretical error for flavour observables, with a 
similar structure to that of BLOCK FOBS, but with the double precision number for the 
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value of the observable replaced by two double precision numbers for the minus and plus 
uncertainties. 

In a similar way, for every block Fname, a corresponding block containing the errors, 
BLOCK FnameERR, can be defined. 

BLOCK FOBSSM 

The block BLOCK FOBSSM contains the SM values of the flavour observables in the same 
format as in BLOCK FOBS. These SM values may be helpful as a reference for comparison. 

BLOCK FDIPOLE 

The block BLOCK FDIPOLE contains the predictions for the electric and magnetic dipole 
moments. The standard for each line in the block should correspond to the FORTRAN 
format 

( lx , 1 10 , 3x , 1 1 , 3x , 1 1 , 3x , IP , E16 . 8 , OP , 3x , ' # ' , lx , A) , 

where the first ten-digit integer should be the PDG code of a particle, the second integer the 
type (electric or magnetic), the next integer the model (SM, NP or SM+NP, as in FWCOEF) 
and finally the last double precision number the value of the moment. The electric dipole 
moments must be given in e.cm unit. 

The PDG codes for the nuclei follows the PDG particle numbering scheme [21] : "Nuclear 
codes are given as 10-digit numbers ±10LZ Z Z AAAI . For a (hyper)nucleus consisting of 
rip protons, n n neutrons and A's, A = n p + n n + tia gives the total baryon number, 
Z = n p the total charge and L = n\ the total number of strange quarks. / gives the isomer 
level, with 1 = corresponding to the ground state (...). To avoid ambiguities, nuclear 
codes should not be applied to a single hadron, like p, n or A , where quark-contents-based 
codes already exist." As an example, the PDG code of the deuteron is 1000010020. 

The types of the moments are defined as follows: 

1 : electric 

2 : magnetic 

The electric and magnetic moments can be given for the SM only, the New Physics (NP) 
only or the total contributions of both SM and NP, using the following definitions: 

: SM 

1 : NP 

2 : SM+NP 

BLOCK FPARAM 

The block BLOCK FPARAM contains process dependent variables, such as form factors, shape 
functions etc., for a specific decay. The decay should be defined as in BLOCK FOBS, but with 
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the type of the observable replaced by a user defined number for the parameter. Here it is 
essential to describe the variable as a comment. 

3 Conclusion 

The interplay between collider and flavour physics is entering a new era with the start-up 
of the LHC In the future more and more programs will need to be interfaced in order to 
exploit maximal information from both collider and flavour data. Towards this end, an 
accord will play a crucial role. The accord presented here specifies a set of conventions 
in ASCII file format for the most commonly investigated flavour-related observables and 
provides a universal framework for interfacing different programs. 

The presented accord will be further developed according to the user feedbacks and 
special needs of specific codes. Some discussions in this direction were held during Les 
Houches 2011 workshop [21], which led to a few clarifications and additions (included in 
the current manuscript). 

As the number of flavour related codes keeps growing, the connection between results 
from flavour physics and high pt physics is becoming more relevant to the disentangling of 
the underlying physics model. Using the lessons learnt from the SLHA, we hope the FLHA 
will prove useful for studies related to flavour physics. 

Acknowledgements 

We thank the organisers of the Les Houches 2009 workshop, where this work was started, for 
their hospitality and support. Thanks also to all those who helped us improving the accord 
by their useful comments, and in particular to G. Belanger, A. Belyaev, F. Boudjema, J.S. 
Lee and A. Pukhov. This work was supported in part by the European Community's Marie 
Curie Research Training Networks under contracts MRTN-CT-2006-035505 'Tools and Pre- 
cision Calculations for Physics Discoveries at Colliders', and MRTN-CT-2006-035606 'MC- 
net'. The work of S.H. was partially supported by CICYT (grant FPA 2007-66387). The 
work of T.G. is supported in part by the Grant-in-Aid for Science Research, Japan Society 
for the Promotion of Science, No. 20244037. 



15 



A The PDG Particle Numbering Scheme 



Listed in the table below are the PDG codes for some important SM mesons. Codes for 
other particles may be found in |21j . 



Name 


PDG code 


Name 


PDG code 




111 


D+ 


411 


7T + 


211 


D° 


421 


p(770)° 


113 


Df 


431 


P(770)+ 


213 


D*+ 


433 


V 


221 


B° 


511 


7/(958) 


331 


B+ 


521 


u(782) 


223 


B*° 


513 


0(1020) 


333 


B *+ 


523 


Kl 


130 


B° s 


531 


K° s 


310 


Bf 


533 


K° 


311 


Bt 


541 


K + 


321 


B *+ 


543 


if* (892) 


313 


J/i?(lS) 


443 


K*+(892) 


323 


T(1S) 


553 


Vc(lS) 


441 


Vb(lS) 


551 



Table 3: PDG codes for most commonly considered mesons. 



B Two-Higgs Doublet Model 

The charged Higgs boson couplings to fermions for the Two-Higgs Doublet Model (THDM) 
can be expressed as 

H+DU : ~ 9 V UD [\ u mu (l - 7 5 ) - \ D m D (l + 7 5 )] , 

2y/2m w ^ 

where U, D and i stand, respectively, for the up-type quarks, the down-type quarks and 
the leptons. The conventions used for the four types of ^-symmetric THDM types, corre- 
sponding to different Yukawa couplings, are given in Table SJ 

The notation and meaning of the different types of model vary in the literature. Some- 
times type Y (III) and type X (IV) are used. In supersymmetry, type III usually refers 
to the general model encountered when the Z 2 symmetry of the tree-level type II model is 
broken by higher order corrections. 
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Type 


x u 


\ D 


x- 


I 


cot /3 


cot /3 


cot /3 


II 


cot/3 


— tan /3 


— tan /3 


III 


cot/3 


— tan /3 


cot /3 


IV 


cot j3 


cot 


— tan /3 



Table 4: Charged Higgs Yukawa coupling coefficients in the ^-symmetric types of the 
THDM. 

C Effective Operators 

Here we give a list of the most relevant effective operators together with their unique two- 
number identifier. 

C.l Effective operators for b — > s transition 

Effective operators relevant to the b — > s transition are 



Oi 

o 2 


= 0(03040405 6161) = 
= 0(03040405 4141) = 
= 0(03050707 4133) = 


(s lfJ ,T a P L c)(crT a P L b) , 
{s^P L c){c^P L b) , 

(s lfl P L b)J2(Q^Q) , 


o 4 


= 0(03050707 6153) = 


(- S1 ,T a P L b)J2(qi"T a q) , 
<? 


o 5 


= 0(03050707 4536) = 


(si^ 2 i m PLb)J2(w^ 2 ^ 3 q) , 




= 0(03050707 6556) = 


(si^l^l^ a PLb)Y,^l^ S T a q) , 
q 


o 7 

Os 


= 0(0305 4422) = -J- 

167T 

= 0(0305 6421) = 9 

167T 


-[sa^(m b P R )b]F^, 
-[sa^(m b P R )T%]G^ . 



C.2 Effective operators for neutral meson mixings 

Effective operators for B® — B® mixing are 

0(01050105 4141) = {d llx P L b){d^P L b) , 
0(01050105 3131) = (dP L b)(dP L b) , 
0(01050105 7171) = (diP^idjPLW) , 
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0(01050105 3132) = (dP L b)(dP R b) , 

0(01050105 7172) = (^P^'X^-P^) , (5) 

and those with opposite chiralities. Operators for B° s — B®, K° — K° and D° — D° mixings 
are defined in the same way. 

C.3 Effective operators for lepton flavour violations 

The effective operators for ALF = 1 lepton flavour violating processes are as follows, 

• n — > e transitions (/x — > e 7, \i — > eee and \i — e conversion in a muonic atom): 



0(1311 4322) 

0(13111111 3131) 
0(13111111 4141) 
0(13111111 4142) 

0(13110101 3131) 
0(13110101 3132) 
0(13110101 4141) 
0(13110101 4142) 
0(13110101 4343) 

0(13110202 3131) 
0(13110202 3132) 
0(13110202 4141) 
0(13110202 4142) 
0(13110202 4343) 

0(13110303 3131) 
0(13110303 3132) 
0(13110303 4141) 
0(13110303 4142) 



m^na^ 'P L e)F flu , 

(jlP L e)(eP L e) , 
(^P L e)(eYP L e) , 
(jl llM P L e)(e^P R e) , 

(fiP L e)(dP L d) , 
(fLP L e)(dP R d) , 
{fl^P L e)(d^P L d) , 
(ft lfl P L e)(drP R d) , 
{jia^P L e){da^P L d) , 

(fiP L e)(uP L u) , 
(jlP L e)(uP R u) , 
(fL-f ll P L e)(u^P L u) , 

{^P^ua^u) , 

(fiP L e)(sP L s) , 
(fxP L e)(sP R s) , 
(fl^P L e)(s^P L s) , 
{m,P L e){s^P R s) , 
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0(13110303 4343) = (pa^ P L e)(sa^P L s) . (6) 

We also define operators with opposite chiralities by replacing Pi and Pr with each 
other in (JBj). 

• r — > fi transitions (r — )■ ^7, r — )■ /i/i/i, r — )■ jie + e~ and r — )■ /ihadrons): 



0(1513 4322) 


= m T (fa^P L fi)F^ , 


0(15131313 3131) 


= (fP L fi)(ftP L fi) , 


0(15131313 4141) 


= (f lfl p L n)(nYP L fi) , 


0(15131313 4142) 


= (f-y (t P L fi)(ji'fP R fi) , 


Dfi ^1 31 1 1 1 31 31 "\ 




^1 31 1 1 1 31 39"\ 




0(15131111 4141) 


= (r 7 , t P L /i)(e 7 ' t P i e) , 


0(15131111 4142) 


= (f^P^ieYPRe) , 


0(15131111 4343) 


= (fa^P Lf i)(ea^P L e) , 


n('\ c;i 3ni m 31 31 

Ly^lOlOUlUl OIOIJ 


— f^r Pr n\( rlP* rl\ 


1 ^1 3m m 31 39^ 


— (rP r lAt^P^rl} 


0(15130101 4141) 


= (f^P Lf i)(dYP L d) , 


0(15130101 4142) 


= {f ltl P L ll){dYP R d) , 


0(15130101 4343) 


= (fa^PL^ida^PLd) , 


0(15130202 3131) 


= (fP L n)(uP L u) , 


0(15130202 3132) 


= (fP L fi)(uP R u) , 


0(15130202 4141) 


= (fj li P L fJ,)(u'y li P L u) , 


0(15130202 4142) 


= (f lfl P L fi)(uYP R u) , 


0(15130202 4343) 


= (fa^PL^iua^PLu) , 


0(15130303 3131) 


= (fP Lf i)(7sP L s) , 


0(15130303 3132) 


= (fP Lf i)(sP R s) , 


0(15130303 4141) 


= (f^P L fi)(s^P L s) , 



19 



0(15130303 4142) = (f^P Lf i)(s^P R s) , 

0(15130303 4343) = (fa^ P L fi)(sa^P L s) . (7) 
Operators with opposite chiralities are also defined. 

• We define operators for r — > e transitions by replacing the muon (code 13) and the 
electron (code 11) fields with each other in (jTJ). 

Effective operators for ALF > 1 leptonic tau decay r + — > \i~ e + e + are 

0(15111311 3131) = (fP L e)(flP L e) , 
0(15111311 4141) = (f llM P L e)(fl^P L e) , 

0(15111311 4142) = {r^P L e)Qi'fP R e) , (8) 

and those with opposite chiralities. Operators for r + — > e~ /i + /i + are defined by replacing 
the muon and electron fields. 

For the processes in which both lepton and quark flavours are violated, such as r — > \x K , 
B° — > /j, e and so on, relevant operators follows, 

0(03011513 3131) = (sP L d)(fP L fx) , 
0(03011513 3132) = (sP L d)(fP R ji) , 
0(03011513 4141) = (s^P L d){f^P LH t) , 
0(03011513 4142) = (s^P L d)(f^P R p) , 

0(03011513 4343) = (sa^ u P L d)(fa^P L ^i) . (9) 

Definitions of operators with opposite chiralities and/or different quark/lepton flavour com- 
binations are straightforward. 

D Decay constants 

The decay constant fp of a pseudoscalar meson P can be defined as: 

(0\qi^Q\P(p)) = -ifpp 11 , (10) 
for q ^ Q quark contents (P = 7r ± , K, D, B). For ir°, rj and rf, we define: 

-L(0\u^ l5 u-dYl 5 d\7T (p)} = -iW, (11) 

-L(0\u^ 5 u + dYl 5 d\^\p)) = -iffaif, (12) 
(0\sn 5 s\ V (>\p)) = -i/;/, (13) 
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assuming isospin symmetry. Other possible choices for 77 and 77' may be: 

-L(0|w 7 At 75« + dri,d-2sn,s\^\p)) = -if^, (14) 

^(0\uY l5 u + dri 5 d + sYl 5 s\^'\p)) = -if^p*. (15) 

In addition, the following matrix elements are defined: 

(m q + m Q )(0\q l5 Q\P(p)) = ih P , (16) 

(m u + m d ) — (0|m7 5 m- dy b d\^{p)) = ih w , (17) 

(m u + m d )-^=(0\uy 5 u + d-f 5 d\r] {,) (p)) = ih 9 (l) , (18) 



V2 



2m s (0|s 7 5^ (,) (p)> = ih'w . (19) 



The parameters h P may be unnecessary except for 77 and 7/ since they can be written in 
terms of other quantities such as h n = rn^f^ etc. h q ^ t) do not satisfy relations of this kind 
due to the contributions of anomaly terms. 

Decay constants of a vector meson V, whose quark content is qQ (such as and K*), 
are defined by the following matrix elements, 

(0\q^Q\V(p)) = m v f v e^, (20) 
(0\q*VQ\V(p)) = ifW^-V^), (21) 

where e M is the polarisation vector of V. f p ,u,<t> i n the "ideal mixing" limit are defined as: 



1 
1 

-.{[)\uYu 

(0\*y"s\<t>(p)) = m+Ue*. (24) 



— (0\uY l u-dYd\p°(p)} = m p f p e», (22) 

M0\u^u + d^d\u(p)) = mj^, (23) 
V2 



fju 4> are a ^ so defined with the same flavour combinations. It is possible to define decay 
constants of u and 4> as 

^=(Q\u^u+d^d\u{cf){p)) = m^f^f^ (25) 

(0\sYsM<P)(p)) = m^f^, (26) 



or 

1 



r _{d\u^u + d^d-2s^s\u{<j>){p)) = m uW f» w e^ (27) 
-^=(0\u^u + d^d + s^ slu^ip)) = m^/^e". (28) 
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E Bag parameters 



We define the bag parameters -81,2,3,4,5 for B d — B d mixing matrix elements as follows, 

9 

* 2 fi 
3 mfl S /B S J 



(B° d \(d^P L b)(drPLb)\B° d ) = ImUloB, , 



(B° d \(dP L b)(dP R b)\B° d ) = (~^-) B 4 , 

z d d \rrib + m d I 



(B^P^djPnb^B ,) = B 5 . (29) 

The renormalisation group invariant bag parameter, Bg, can also be defined in function of 
B\. Bag parameters for B° s — B®, K° — K° and D° — D° mixings are defined in the same 
way. 



F Meson mixings 

Meson mixing is assigned observable type 7 in BLOCK FOBS. In the FLHA, we assume the 
standard definition for the meson mixings. The oscillation frequency of Q° q and Q Q q mixing 
is characterised by the mass difference of the heavy and light mass eigenstates [25] : 

AM q = M| - Ml = 2Re ^/(Mf 2 - ^(Mf* - ^), (30) 

where M12 and Ti2 are the transition matrix elements from virtual and physical intermediate 
states respectively. For the kaon systems this gives: 

AM K = 2ReM 12 . (31) 

In the case of AB = 2 transitions since | 1 "C IM12I, one can write: 

(B q >\U* B = 2 \B° q } = 2M B M? 2 , (32) 

where Ms q is the mass of B q meson and 

AM q = 2\M q 12 \. (33) 

The quantity given in the block is AM q and the units are fixed to 1/ps. We use the PDG 
number of the oscillating mesons for the parent and the flipped sign for the daughter^ The 
number of daughters is fixed to 1. For example, the B s — B s oscillation frequency is given 
as: 

2 Exceptions to this rule are possible, for instance, in the case of K-K mixing, where the states K\ ong 
and -Kshort with their PDG numbers 130 and 310, respectively, can be used. 
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Block FOBS # Flavour observables 

# ParentPDG type value q NDA ID1 ID2 ID3 . . . comment 

531 7 1.9e01 1 -531 # Delta M_s 

Similarly, the corresponding SM values and the errors can be given in BLOCK FOBSSM 
and BLOCK FOBSERR, respectively. Note that the matrix elements (which are not physical 
observables) cannot be given in this block. Such quantities can be expressed in terms of 
Wilson coefficients, decay constants, bag parameters, etc. which are defined in the FLHA, 
or they can be given in a user defined block. 



G FLHA sample file 

The following is an example of a generic FLHA file. For simplicity, we have chosen a SUSY 
model, but instead of giving the complete SLHA spectrum, we here present it as a generic 
FLHA model. Therefore we use FMODSEL to specify the model switches rather than using 
MODSEL, which would necessitate repeating the complete SLHA spectrum. For examples of 
the latter, see the SLHA manuals [18j[T9]. Some lines in the blocks FOBS and FOBSMS are 
broken for better readability. 



# Superlso output in Flavour Les Houches Accord format 
Block FCINFO # Program information 

1 SUPERISD # flavour calculator 

2 3.2 # version number 
Block FMODSEL # Model selection 

1 1 # Minimal supergravity (mSUGRA,CMSSM) model 



6 


1 


# Quark flavour 


is violated 


Block SMINPUTS # Standard 


Model inputs 


1 


1 


,27934000e+02 


# 


alpha_em~ (-1) 


2 


1 


. 16637000e-05 


# 


G_Fermi 


3 


1 


. 18400000e-01 


# 


alpha_s(M_Z) 


4 


9 


. 11876000e+01 


# 


m_{Z}(pole) 


5 


4 


. 19000000e+00 


# 


m_{b}(m_{b» 


6 


1 


,72900000e+02 


# 


m_{top}(pole) 


7 


1 


,77700000e+00 


# 


m_{tau}(pole) 


Block MINPAR 


# SUSY breaking input parame 


1 


8 


.00000000E+02 


# 


m_0 


2 


5 


.00000000E+02 


# 


m_{l/2> 


3 


2 


.00000000E+01 


# 


tan(beta) 


4 


1 


.00000000E+00 


# 


sign(mu) 


5 


-5 


.00000000E+02 


# 


A_0 



Block FMASS # Mass spectrum in GeV 

#PDG_code mass scheme scale particle 

3 1.04000000e-01 1 2 . 00000000e+00 # s (MSbar) 

5 4.68474767e+00 3 # b (IS) 

211 1.39600000e-01 # pi+ 
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313 
321 
421 
431 
521 
531 
Block FLIFE 



8.91700000e-01 
4.93700000e-01 
1.86960000e+00 
1.96847000e+00 
5.27917000e+00 
5.36630000e+00 
# Lifetime in 



# K* 

# K+ 

# DO 

# D_s+ 

# B+ 

# B_s 

sec 
particle 

# pi+ 

# K+ 

# D_s+ 

# B+ 

# B_s 

Block FCONST # Decay constant in GeV 

#PDG_code number decay_ const ant scheme scale particle 

431 1 2.48000000e-01 # D_s+ 

521 1 1.92800000e-01 # B+ 

531 1 2.38800000e-01 # B_s 

Block FC0NSTRATI0 # Ratio of decay constants 
#PDG_codel code2 nbl nb2 ratio scheme scale comment 

321 211 1 1 1. 19300000e+00 # f_K/f_pi 
Block FPARAM # Process dependent parameters 



#PDG_code 
211 
321 
431 
521 
531 



lifetime 

2.60330000e-08 

1.23800000e-08 

5.00000000e-13 

1.63800000e-12 

1.42500000e-12 



ParentPDG 


number value 


NDA 


ID1 


ID2 


ID3 


5 


1 


5.80000000e-01 


2 


3 


22 




521 


1 


4.60000000e-01 


3 


421 


-15 


16 


521 


2 


1.03000000e+00 


3 


421 


-15 


16 


521 


3 


1.17000000e+00 


3 


421 


-15 


16 


521 


1 


2.68000000e-01 


2 


313 


22 





Block FWC0EF Q= 1 . 60846e+02 



. comment 

# C in b->s gamma 

# Delta(w) in B+->D0 tau nu 

# G(l) in B+->D0 tau nu 

# rho"2 in B+->D0 tau nu 

# Tl in B->K* gamma 
# Wilson coefficients at scale Q 



id 




order 


M 


value 


comment 


03040405 


6161 


00 


2 


0. 


00000000e+00 


# 


CI" 





03040405 


4141 


00 


2 


1 


,00000000e+00 


# 


C2" 





0305 


4422 


00 


2 


-1. 


,53496321e-01 


# 


C7" 





0305 


6421 


00 


2 


-9, 


,51462419e-02 


# 


C8" 





03040405 


6161 


01 


2 


2. 


,33177662e+01 


# 


ci- 


1 


03040405 


4141 


01 


2 





,00000000e+00 


# 


C2" 


1 


03050707 


6153 


01 


2 


5 


29858390e-01 


# 


C4" 


1 


0305 


4422 


01 


2 


-4, 


27127914e-01 


# 


C7" 


1 


0305 


6421 


01 


2 


-1. 


,06024128e+00 


# 


C8" 


1 


03040405 


6161 


02 


2 


3. 


,08457152e+02 


# 


CI" 


2 


03040405 


4141 


02 


2 


4 


,91587899e+01 


# 


C2" 


2 


03050707 


4133 


02 


2 


-6, 


,99555993e+00 


# 


C3" 


2 


03050707 


6153 


02 


2 


1. 


,25585165e+01 


# 


C4" 


2 


03050707 


4536 


02 


2 


8 


,73907482e-01 


# 


C5" 


2 


03050707 


6556 


02 


2 


1 


63857653e+00 


# 


C6" 


2 


0305 


4422 


02 


2 


8. 


,40063580e-01 


# 


C7" 


2 


0305 


6421 


02 


2 


-4, 


,61747635e+00 


# 


C8" 


2 
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Block FDIPOLE # Electric and Magnetic dipole moments 
# PDG_code type M value comment 

13 2 1 5.00716589e-10 

Block FOBS # Flavour observables 



ParentPDG 


type value 




q 


5 


1 


2 


83249483e 


-04 


o 


521 


4 


8 


78354299e 


-02 





531 


1 


3 


44359657e 


-09 





521 


1 


9 


89829667e 


-05 





521 


2 


9 


81528344e 


-01 





431 


1 


5 


10125129e 


-02 





431 


1 


5 


23481907e 


-03 





521 


1 


6 


71755629e 


-03 





521 


11 


2 


96312243e 


-01 





321 


11 


6 


34141764e 


-01 





321 


12 


9 


99918554e 


-01 
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# BR(b->s gamma) 
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BR(B->D0 e nu) 

# BR(K->mu nu)/ 
BR(pi->mu nu) 

# R_123 

ID3 . . . comment 

# BR(b->s gamma) 

# Delta0(B->K* gamma 

# BR(B_s->mu+ mu-) 

# BR(B_u->tau nu) 

# R(B_u->tau nu) 

# BR(D_s->tau nu) 

# BR(D_s->mu nu) 
16 # BR(B->D0 tau nu) 
16 # BR(B->D0 tau nu)/ 

BR(B->D0 e nu) 

# BR(K->mu nu)/ 
BR(pi->mu nu) 

# R_123 



References 



[1] F. Mahmoudi, Comput. Phys. Commun. 178 (2008) 745 [arXiv:0710.2067] : Comput. 
Phys. Commun. 180 (2009) 1579 [arXiv:0808.3l44j : Comput. Phys. Commun. 180 
(2009) 1718. 



Code website: http : / /superiso . in2p3 . f r 



[2] G. Degrassi, P. Gambino and P. Slavich, Comput. Phys. Commun. 179 (2008) 759 
[arXiv:0712.3265j . 



Code website: http : / /slavich . web . cern . ch/ slavich/susybsg 



25 



[3] 
[4] 

[5 

[6 

[7 



[9 

[10 

[11 
[12 

[13 

[14 
[15 



P. Paradisi, talk given at "Interplay of Collider and Flavour Physics, 3rd general meet- 
ing", CERN, December 2009. 

J. Rosiek, P. Chankowski, A. Dedes, S. Jager and P. Tanedo, Comput. Phys. Commun. 
181 (2010) 2180 jarXiv: 1003.4260] . 



Code website: http://www.fuw.edu.pl/susy_flavor 



G. Belanger, F. Boudjema, A. Pukhov and A. Semenov, Comput. Phys. Commun. 180 
(2009) 747 [arXiv:0803.2360j . 



Code website: http://wwwlapp.in2p3.fr/lapth/micromegas 



A. Arbey and F. Mahmoudi, Comput. Phys. Commun. 181 (2010) 1277 
[arXiv:0906.0369j . 



Code website: http://superiso.in2p3.fr/relic 



S. Heinemeyer, W. Hollik and G. Weiglein, Comput. Phys. Commun. 124 (2000) 76 
|hep-ph/9812320|; Eur. Phys. J. C 9 (1999) 343 |hep-ph/9812472] . 
Code website: |http://www. f eynhiggs . de , 

J. S. Lee, A. Pilaftsis, M. S. Carena, S. Y. Choi, M. Drees, J. R. E llis and C. E. M. Wag- 
ner, Comput. Phys. Commun. 156 (2004) 283 |hep-ph/0307377|. 
Code website: |http : //www . hep . man . ac . uk/u/ j slee/ CPsuperH . html 



U. Ellwanger and C. Hugonie, Comput. Phys. Commun. 175 (2006) 290 
lhep-ph/0508022). 



Code website: http://www.th.u-psud.fr/NMHDECAY/nmssmtools.html 



F. E. Paige, S. D. Protopopescu, H. Baer and X. Tata, hep-ph/03 12045 



Code website: http : / / www . nhn . ou . edu/~ isaj etj 

W. Porod, Comput. Phys. Commun. 153 (2003) 275 [hep-ph/0301101 



Code website: http : //www . physik . uni-wuerzburg . de/~ porod/ SPheno . html 



A. Djouadi, J. L. Kneur and G. Moultaka, Comput. Phys. Commun. 176 (2007) 426 
[hep-ph/02 11331). 



Code website: http : //www . lpta . univ-montp2 . f r/users/kneur/Suspect 



B. C. Allanach, Comput. Phys. Commun. 143 (2002) 30 5jhep-ph/0104145| . 
Code website: http : //pro j ects . hepf or ge . org/ sof tsusy| 



R. Lafaye, T. Plehn and D. Zerwas, hep-ph/0404282 



P. Bechtle, K. Desch and P. Wienemann, Comput. Phys. Commun. 174 (2006) 47 
[hep-ph/0412012 



Code website: |http : //www-f lc . desy . de/ f ittino 



[16] R. Ruiz de Au stri, R. Trotta and L. Roszkowski, JHEP 0605 (2006) 002 
[hep-ph/0602028 



Code website: |http://www.ft .uam.es/personal/rruiz/superbayes 



26 



[17] S. Heinemeyer, talk given at "Interplay of Collider and Flavour Physics, 2nd general 
meeting" , CERN, March 2009. 



Code website: http://cern.ch/mastercode 



[18] P. Skands et al, JHEP 0407 (2004) 036 [hep-ph/ 031 1123] . 

[19] B. Allanach et al, Comput. Phys. Commun. 180 (2009) 8 |arXiv:0801.0045] . 
Website: |http : / /home . f nal . gov/- skands/ slha[ 

[20] T. Hahn, |hep-ph/0408'283l Comput. Phys. Commun. 180 (2009) 1681 



lhep^ph/0605049 



Website: |http : //www . f eynarts . d e/ slha 



[21] K. Nakamura et al. [Particle Data Group Collaboration], J. Phys. G G 37 (2010) 
075021. 

[22] W. Siegel, Phys. Lett. B 84 (1979) 193. 

[23] I. Jack et ai, Phys. Rev. D 50 (1994) 5481 |hep-ph/940729T] . 



[24] |http : //phystev . in2p3 . f r/Houches201 i 



[25] A. J. Buras, hep-ph/9806471 



27 



